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Abstract 

We study the evolution of abelian-Higgs string networks in numerical 
simulations. These are compared against a modified velocity-dependent 
one scale model for cosmic string network evolution. This incorporates 
the contributions of loop production, massive radiation and friction to 
the energy loss processes that are required for scaling evolution. We find 
that the loop distribution statistics in the simulations are consistent with 
the long-time scaling of the network being dominated by loop production. 
For an oscillating sinusoidal perturbation, we also demonstrate that the 
power emitted into massive radiation decays strongly with wavelength. 
Putting these observations together and extrapolating, we believe there is 
insufficient evidence to reject the the standard picture of string network 
evolution in favour of one where direct massive radiation is the dominant 
decay mechanism, a proposal which has attracted much recent interest. 


1 Introduction 

Vortex-string networks are important in a variety of contexts, whether in con¬ 
densed matter physics or cosmology [0. If we are to obtain a quantitative de¬ 
scription of these networks, then we must also properly understand their ‘scaling’ 
evolution as well as the decay mechanisms which maintain it. The abelian-Higgs 
model, a relativistic version of the Ginzburg-Landau theory of superconductors, 
provides a convenient testbed for developing detailed models for this evolution. 
On the one hand, the relatively simple field theory can be studied directly in 
three-dimensional simulations. On the other, a straightforward reduction to a 
one-dimensional effective theory—the Nambu action—can also be studied nu¬ 
merically, though over a much wider dynamic range. 

In a cosmological context, a rather simple ‘one-scale’ model of string evo¬ 
lution has emerged ||. || which appears to successfully describe the large-scale 
features of an evolving string network Q, ^|, though with subtleties remaining on 
smaller scales. In this simple model, the average number of long strings in a hori¬ 
zon volume remains fixed as it expands, a rapid dilution made possible through 
reconnections resulting in loop production. The loops, in this standard picture, 
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oscillate relativistically and decay through gravitational radiation. The subtlety 
here concerns small loop creation; gravitational radiation backreaction effects 
should act on the long string network eliminating small wavelength modes, thus 
setting a minimum loop creation size||, This backreaction length ‘scales’ 
with the horizon size (for GUT-scale strings it should be approximately 10” 4 f), 
and so loop sizes should also be scale-invariant, albeit tiny and, as yet, not 
adequately probed by Nambu string simulations. 

Recently this standard picture for network evolution has been questioned on 
the basis of abelian-Higgs field theory simulations!^. The authors suggest that 
the primary energy loss mechanism by long strings is direct massive radiation, 
rather than loop creation. This is contrary to qualitative expectations that 
the presence of a large mass threshold should exponentially suppress massive 
particle production for any long wavelength oscillatory string modes, that is, 
those much larger than the string widthQ. Evidence in support of this claim 
rests primarily on the study of large amplitude oscillations of a single string and 
network simulations in which the loop density is observed to be low. 

The aim of the present work is to consider these issues by taking a more 
detailed look at radiation from an oscillating string and the available decay 
mechanisms for an evolving network. In high resolution and low noise simu¬ 
lations of a perturbed string, we are able first to demonstrate that qualitative 
expectations for massive radiation suppression are correct. Next, by analyti¬ 
cally modelling the small-scale results from field theory network simulations, 
we are able to argue that these can be sensibly and consistently extrapolated 
to the large lengthscales relevant for Nambu simulations and cosmology. Al¬ 
though complex nonlinear processes are at work on small-scales in the field 
theory simulations, loop and ‘protoloop’ production already appears to be the 
be the dominant energy decay mechanism. 


2 Massive radiation from an oscillating string 

As stated above, we shall be considering strings in the Abelian Higgs model, 
the simplest producing gauged vortex-lines. The Lagrangian density is 

C = (D^riD^) - -F^ - - T7 2 ) 2 . (1) 

The covariant derivative D M acts on <j> as d^ — ieA^. The Higgs field mass is VA 77 
and that of the vector field is \f2er\. After rescaling the only free parameter in 
the model is the ratio of these two masses, which we take to be unity (specifically, 
A/2 = r\ = e = 1). This corresponds to the Bogomol’nyi limit for vortices in two 
dimensions. The system of evolution equations resulting from (([!])) is solved 
using the standard lattice gauge theory methods proposed by Myers et al. [||| 
(we also always adopt the temporal gauge A 0 = 0). Initial data is evolved 
forward in time using a leap-frog discretization of Hamilton’s equations with a 
short time-step. Throughout our simulations the boundary conditions in use 
are periodic up to a gauge transformation in all cases and exactly periodic for 
network simulations. 

First, we consider a perturbed string lying in the ^-direction and oscillating 
in a periodic box of sidelength L. The starting configuration for each simulation 
is based on a simple ansatz for which the core position varies sinusoidally along 
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the string length with amplitude y = A sin(27rz/L). We use profile functions 
obtained by solving the one dimensional field equations for a cylindrical string 
numerically. If, as we suggest, the perturbed string is weakly radiating for large 
L then it is necessary to distinguish the ’true’ radiation from that due to the 
spurious modes introduced by inappropriate initial conditions. We achieve this 
by relaxing the gauge links in our lattice system to their energetic minimum, 
while fixing the Higgs field, and hence the position of the string core (this 
removes the worst modes). Further spurious modes associatied with the Higgs 
field modulus are then removed by a releasing all the fields for a short period 
to evolve via ‘gradient flow’, that is, purely first-order dissipative evolution. 



Length 


Figure 1: The power per unit length emitted by a oscillating strings as a 
function of length L, initially given a fixed relative amplitude e = 2i tA/L = 0.9. 
The fall-off with increasing length is not consistent with the L -1 required for 
scaling. The dashed line is a best least-squares fit to exponential fall-off. 

To calculate radiative power we look at the rate of change of energy outside 
a region of fixed radius around the position of the unperturbed string. We 
take the average initial energy increase as the radiation rate and the effect of 
varying the string length is illustrated in Figure Here, the relative amplitude 
for the initial data was kept at a large fixed value, e = 27 tA/L = 0.9; it is 
important that this value is below unity because, for e > 1.0, the perturbed 
strings correspond to degenerate relativistic Nambu configurations which radiate 
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pathologically An exponential trend as a function of length L is apparent, 
which can be approximately fitted by P oc exp (—L/L m ) with L m « 4.8. For 
smaller amplitude, there is an even stronger dependence on L, so we can be 
confident that L m < 4.8 for e < 0.9 (as we describe in more detail elsewhere 
[fTT]]). Certainly the decay of emitted power between the shortest and longest 
L shown here is consistent only with a higher power law than is required for 
scaling, that is, P oc L~ n with n » 1. For L significantly larger than shown in 
Figure [I], the minute power loss becomes indistinguishable from the background 
numerical noise inherent in the simulations. Note that the actual spectral nature 
of this radiation and its relationship to the string perturbation lengthscale L is 
discussed at length elsewhere Jn) . 

The strong scale-dependence of massive radiation we observe is not consis¬ 
tent with the behaviour found in ref. JtJ , where they suggest that the power per 
unit length scales merely as P oc L~ x . We can see two possible explanations 
for this clear discrepancy: First, the results reported in ref. j7]J are for only a 
single relative perturbation amplitude e = ty. This is a strongly nonlinear and 
degenerate regime for which pathologically strong radiation is expected. Sec¬ 
ondly, for large amplitudes the unrelaxed initial ansatz in ref. [Q is inaccurate, 
particularly for the gauge fields. As we have observed in our own simulations, 
this can lead to spurious modes from the initial relaxation which swamp the 
radiation due to the string oscillations. 

3 String network evolution 

We have performed network simulations using periodic cubic lattices of side 
length 250 and greater. The principal physical parameter in each simulation is 
the initial correlation length. To establish a suitable network configuration we 
take an initial configuration with a flat initial power spectrum of fluctuations 
in the Higgs field, centered around zero, and zero gauge field. This is evolved 
forward in time with Hamiltonian dynamics to increase the correlation length 
and then acted upon by a short period of gradient flow evolution to reduce the 
initial inter-string energy density to a negligible value, since the one-scale model 
requires that this be zero initially. Subsequently the simulation is evolved using 
Hamiltonian dynamics. The Gauss constraint is satisfied initially as we start 
from rest, and is preserved by the equations of motion. 

The time-step is chosen to be sufficiently small to allow the Hamiltonian to 
be conserved within 1 per cent over the course of a run. The choice of spatial 
lattice spacing is constrained by the desire to meet two conflicting criteria: (i) 
to simulate a volume which is orders of magnitude larger than the string width, 
and (ii) to accurately represent the continuum theory. Throughout we have 
chosen a lattice size of 0.5, which is as large as we feel is reasonable given that 
at larger spacings there is a significant potential barrier associated with the 
lattice and that for oscillating strings of length « 15 we have observed greatly 
increased radiation using a larger lattice spacing. 

1 For a relative amplitude e > 1. a sinusoidal Nambu string develops ‘lumps’, that is, finite 
regions of string which pile up at a point moving at the speed of light. Not surprisingly, 
radiative backreaction in a field theory simulation is severe in such regions, leading to the loss 
of a constant length of string in the first oscillation; initially, then, we would have a power 
per unit length P ~ L~ l . We do not expect such perfectly degenerate string configurations 
in a cosmological context. 
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To characterise the network configuration at a given point in time we assign 
positions of zeros of the Higgs field to lattice plaquettes according to the winding 
of the Higgs phase around each plaquette. This allows us to calculate correlation 
lengths and loop distribution statistics ||. When combined with similar data 
from a nearby time-step we can also estimate the velocity of each string segment. 
The procedure for calculating velocities is relatively vulnerable to numerical 
errors due to uncertainty in the ‘true’ position of the string network and the 
sensitivity of measured velocities to this. 

In Figure]^ we see general consistency with the results of Figure 1 of Vincent 
et al. . In both cases the network correlation length grows approximately linearly 
throughout the simulations. The most notable deviations from linearity is seen 
in the run at largest lattice spacings of ref. [Q, where the rate of growth of the 
correlation length with time is ~ 0.5, some 50 per cent greater than in the rest 
of the simulations. We speculate that this increased rate of growth is caused 
by lattice effects. The initial brief burst of growth in correlation length seen 
by Vincent et al. in several simulations is not seen in our simulations and may 
arise from the slightly different procedures used to set up the initial network 
configurations. 



Figure 2: The spatial correlation length L as a function of time for a series of 
simulations with a grid size of 250 3 . A velocity dependent ‘one-scale’ model is 
used to fit the data (solid lines). 
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4 Analytic modelling of network evolution 

If we regard the one-dimensional Nambu action as providing a satisfactory first 
approximation to string evolution then we can employ a velocity-dependent 
‘one-scale’ model to describe the evolution of a string network g]. The long 
string network energy density p x is susceptible to three possible energy loss 
mechanisms, friction, loop production and direct massive radiation which are 
phenomenologically summarized in the following averaged evolution equation 

@ : 

dpoo V Poo CVPoo / i f ^ 

-dT = — + — + p ( 2 ) 

where Lf is the friction length, c is the loop chopping efficiency, and /(L/L m ) 
is proportional to the power per unit length of massive radiation. Employing 
the definition for the correlation length p = p/L 2 and noting a supplementary 
equation for the rms velocity v we obtain ||] 
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We provide some explanation for the origin of these terms below, but note 
that the first term in the velocity equation is due to the acceleration of strings 
with a typical curvature radius L (with v l/\/2 for free relativistic motion). 
Given the difficulties in measuring and normalizing this velocity in the field 
theory simulations, we also allow for lattice discretization effects through the 
phenomenological parameter a; this term takes a form which could also incor¬ 
porate momentum losses due to loop production. 

A potentially important energy loss mechanism for the string network is 
friction through interactions with background energy density p. For the simu¬ 
lations described here with dissipative initial conditions, this background arises 
dynamically through the decay of the string network into massive particles, that 
is, p = Poo,i — Poo ■ Characterising this in terms of a friction lengthscale above 
which friction dominates, we have Lf = p(3/p = (3(l/L 2 — 1/L 2 ) -1 (for L — > oo 
in flat space, Lf — > f3L 2 . We can bound the friction coefficient (3 by beginning 
with a high density string network and assuming a late-time friction-dominated 
regime. An example is the lowest curve in Figure ||, where friction affects the 
scaling behaviour, most obviously through the velocity in Figure |j; this pro¬ 
vides a lower limit of (3 > 0.2. From this we can infer that friction provides less 
than 10% of the energy losses for the duration of the simulations in Figure |], 
if the initial correlation length satisfies L\ > 20. Hence, assuming friction to be 
insignificant for the lowest density simulations, we can make a more accurate 
estimate /3 ~ 0.3. 

Loop production is caused by the intersections and self-intersections of long 
strings. For the simulations least affected by friction, let us assume the standard 
picture with energy losses predominantly due to this loop production. In this 
case, from a fit to the correlation length and velocities we can obtain the asymp¬ 
totic value c ss 0.75 (with some small momentum losses a ss 0.45). This value of 
c is consistent with previous determinations of the loop chopping efficiency in 
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flat space simulations |l 2 | ; note that the velocity-dependence in our definition 
implies that their c = c/v). 

Direct massive radiation will provide some energy losses in the early evo¬ 
lution, but this should diminish as the correlation length increases, consistent 
with the behaviour of the oscillating string discussed in section 2. From Figure 
we can see that there is considerable evidence for an additional scale-dependent 
initial energy loss mechanism, since the initial slopes for the smaller correlation 
lengths are clearly steeper than those with L{ > 15. From a fit to Figure ||, 
the function describing the energy losses into massive radiation should behave 
approximately as f(L/L m ) ss 7 exp (—L/L m ) with L m ~ 4.8 for the present pa¬ 
rameters. The overall strength 7 can only be rougly estimated from the initial 
slopes in Figure[2]to yield 7 ~ 0.15. This implies, consistent with Figure that 
massive radiation provides a strong initial contribution only for L < 10 and, 
subsequently, it is rapidly curtailed. 

What we obtain finally is the self-consistent fit of our analytic model to the 
correlation length and velocity evolution shown in Figure |j. Asymptotically this 
is quantitatively in agreement with the standard picture of long string evolution, 
and it also reproduces predicted qualitative features. However, if we choose to 
ignore the loop contribution as in ref. |t| and we assume that massive radiation 
losses are scale-invariant (/ oc A -1 ), then it is also possible to fit the asymptotic 
data with this analytic model, although the initial qualitative features would not 
be as consistent. The question of whether the dominant loop loss interpretation 
is correct, therefore, must hinge on more precise information about the loop 
distribution and how it evolves in time. 


5 Competing energy loss mechanisms 

In using these simulations to determine the dominant decay mechanisms for a 
cosmic string network we face a rather obvious difficulty. If the standard cosmo¬ 
logical picture of loop production is correct, then we know from high resolution 
Nambu string simulations that the typical loop creation scale £ relative to the 
horizon is a = £/f < 10 “ 3 , while radiation backreaction estimates suggest that 
a ~ 10 -4 for GUT-scale strings. The present field theory simulations, however, 
have a dynamic range which is at least two orders of magnitude poorer, so we 
cannot realistically hope to probe these small-scale regimes. In consequence, 
we might be surprised to be able to identify any loops at all and we certainly 
cannot expect their average creation size £ to ‘scale’ relative to L. As we ob¬ 
serve numerically, most loops which are created have radii comparable to the 
string thickness, with £ ~ 7 r almost constant throughout. Indeed, simulation 
visualisations indicate that most energy is lost via ‘protoloops’, that is, small- 
scale highly nonlinear, but coherent, regions of energy density (as illustrated in 
Figure [|). Consistent with the Nambu simulations, this ‘protoloop’ production 
occurs—like small loop formation—in high curvature regions where the strings 
collapse and become very convoluted [pdf . 

A proportion of these ‘protoloops’ have sufficient topology to be identified as 
loops by our numerical diagnostics, so we can estimate the energy loss via this 
pathway relative to other mechanisms. Because of dynamic range limitations, all 
small loop trajectories in these simulations are self-intersecting and loops decay 
almost as fast as allowed by causality. Their time-dependent length, therefore, 



Figure 4: Four stages of the formation of ‘protoloops’ in the field theory 
simulations. : A highly curved region of string (i) collapses to form nonlinear 
‘lumps’ in the energy density (ii) which are lost (iii,iv) when the the string’s 
topological winding intercommutes or annihilates in this region. Energy density 
contours are plotted. 
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can be given by £ « £ c — 4 (t — t c ) in the range 0 < t — t c < £ c /4 with the loop 
creation length £ c and time t c . Hence, at any one time, the string energy loss 
through loop decay can be simply approximated by pn ~ —4 pri£, where ng is 
the loop number density.^] For linear scaling to pertain, a dominant energy loss 
mechanism must behave as p oc t ~ 2 3 , implying n oc t~ 3 so the loop energy density 
would behave as pi = p£ne. oc Poo/t (with £ constant). Thus, contrary to one 
of the conclusions in ref. (t| , loops can be an important decay mechanism even 
though their relative contribution to the total string energy density is falling 
rapidly as t~ l . 

Figure illustrates the relative small loop contribution to the overall energy 
density losses throughout the particular simulation beginning with L[ = 15. 
This is compared with the estimated analytic fit for friction, loop losses and 
massive radiation. We can observe that the measured loop energy losses grow 
steadily towards the analytic loop contribution, which for these simulations 
we assume must include ‘protoloops’ both with topology and without it. We 
can see that the proportion of loops—the topological ‘protoloops’—gradually 
grows to meet or even overtake the analytic loop contribution by the end of the 
simulation. In fact our analytic calculation deserves closer quantitative scrutiny 
because it is a significant overestimate; the histogram plotted in Figure |H] only 
gives p,ri(L 3 rather than four times this quantity. 

With this clear qualitative trend evident, it is not unreasonable to conjecture 
that an extrapolation by another twenty orders of magnitude to cosmological 
scales will imply that the loop contribution will be completely dominant. As 
the typical string perturbation lengthscale grows and is affected by radiative 
backreaction, it is again reasonable to suppose that the typical loop creation 
size will also grow, becoming many orders of magnitude larger than the string 
thickness. We conclude that these field theory simulations, once we account 
for their small dynamic range, are consistent with the standard picture of long 
string network evolution via small loop production. At the very least, the simu¬ 
lations do not provide compelling evidence that cosmological strings will decay 
primarily through the direct radiation of ultra-massive particles. 


6 Conclusion 

We conclude from these numerical results and their analytic interpretation that 
the standard picture provides a more coherent and adequate model for string 
network evolution than the more radical alternative based on direct massive 
radiationjTj]. However, this is not to suggest that we have provided a complete 
or detailed description of the complex nonlinear processes that underlie network 
evolution on these small scales. Rather, first, we have shown qualitatively that 
loop production is important for network evolution and should become more so 
when extrapolated to large scales. Secondly, we have also demonstrated that 
massive radiation is strongly suppressed for long wavelength modes, implying 
that it is an inadequate decay channel for maintaining network ‘scaling’. We 
are currently investigating both these aspects in more quantitative detail[pT[. 

These results have significant cosmological implications. Our expectation is 
that massive particles will only be produced infrequently in highly nonlinear 

2 Note that this is quite different to the cosmological Nambu limit where loops can be 

non-intersecting and long-lived. 
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Figure 5: The relative contributions of the different components to the overall 
energy density losses of the string network, as predicted by the analytic model 
for the simulation in Figure with L\ = 15. Superposed onto this figure is 
a histogram of energy loss into loop production estimated directly from mea¬ 
surements of the time-varying loop density in this particular simulation. It is 
apparent that the loop energy loss contribution is becoming more important as 
the simulation progresses. 
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string regions, such as at cusps and reconnections. The ensuing flux of cosmic 
rays should be relatively low|l4|, |Ii|. The more recent estimates of the ensuing 
cosmic ray flux from direct massive radiation from strings appear to be overly 
optimistic [US). This work points to the need for caution in making cosmolog¬ 
ical extrapolations from small-scale numerical simulations and to the need for 
further progress understanding string radiation backreaction. 
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